SKIN DEPTH COMPENSATION IN UNDERGROUND BORING APPLICATIONS 



BACKGROUND OF THE INVENTION 

The present invention is related generally to the field of locating and/or guiding an underground boring tool 
using a locating signal which is transmitted through the ground and, more particularly, to a method and associated 
apparatus for locating and/or guiding the boring tool in a way which compensates for skin effect that potentially 
introduces error in locating and/or guiding the boring tool as a result of conductivity of the earth through which the 
locating signal passes. 

Referring to Figure 1, boring tools are typically guided or located by transmitting a dipole field from a dipole 
transmitter which is positioned within the drill head of the boring tool. The locating/dipole field is an oscillating signal 
that is generally emitted from a dipole antenna oriented along the rotational axis of the drill head. Figure 1 illustrates a 
coordinate system including x, y and z axes with a dipole transmitter D at its origin. For a point p, at a radius r from the 
origin, the dipole equations are given as: 
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Where B x , B y and B z represent orthogonal components of the dipole field at point p. The dipole equations are 
recited herein for the benefit of the reader since these equations form a fundamental basis for the use of a dipole field in 
locating applications. One such locating system is described, for example, in US Patent no. 5,337,002 which is 
commonly assigned with the present application. Traditionally, boring tool systems have not used compensation for 
conductivity of the soil even though this conductivity introduces a phenomenon commonly referred to as skin effect. 
While skin effect can result in significant locating errors, applicants submit that prior art systems have not provided such 
compensation, at least in part, since it is perceived in the art that compensation for skin effect is an extremely complex 
proposition. 

What prior art system designers have generally done is to altogether ignore skin effect. This is tantamount to an 
assumption of a non-conducting earth. Accordingly, the electromagnetic field emitted by the magnetic dipole of a 
transmitter into a non-conducting medium (such as air) is described mathematically by the well known cubic law of a 
magnetic dipole (see Figure 1). Unfortunately, however, as a direct result of skin depth, drilling in the earth can produce 
significant deviations from the cubic law when a typical oscillating magnetic dipole field is used. The latter term 
describes a magnetic dipole having a signal strength that varies sinusoidally with time. 
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The present invention provides a highly advantageous and heretofore unseen method and associated apparatus 
which provide compensation for skin effect in underground boring tool applications. 



SUMMARY OF THE INVENTION 

As will be described in more detail hereinafter, there are disclosed herein arrangements, apparatus and 
5 associated methods for skin depth compensation in underground boring applications. Accordingly, in an overall method 
of operating a system in which a boring tool is moved through the ground in a region which includes an electrical 
conductivity characteristic and where the system includes an above ground arrangement for tracking the position of 
and/or guiding the boring tool as the boring tool moves through the ground and in which the system is configured for 
transmitting a locating signal between the boring tool and the arrangement in the region, the improvement comprises 
1 0 compensating for skin depth error by measuring the locating signal such that measurements of the locating signal include 
skin depth error introduced as a result of the electrical conductivity characteristic and, thereafter, using the measurements 
in a way which determines a skin depth corrected position of the boring tool. 

In one aspect of the invention a multi-frequency approach is provided which utilizes measured intensities of the 
O locating field at two or more frequencies to extrapolate a zero frequency value of locating signal intensity. The zero 
frequency value of intensity is then used in position determination. The multi-frequency approach may be used in 
m conjunction with walk-over type locators or with one or more above ground receivers designed for receiving the locating 
-p signal at fixed position(s). In one feature, the multi-frequency approach of the present invention does not require 
H knowledge of earth properties or ground surface geometry. The components of the measured magnetic field intensities 
m of the locating field measured at their selected frequencies contain property and geometry effects and pass them on to 
-20 extrapolated zero frequency values. 

= y ; In another aspect of the invention, certain intensity measurements of the locating signal are used to determine a 

jr value for skin depth to be used during subsequent drilling, these certain measurements being obtained in a calibration 
O procedure by transmitting the locating signal from the boring tool on the surface of the ground to the above ground 
arrangement prior to drilling. 

25 In still another aspect of the invention, a determined value of skin depth is used in one locating scenario with a 

walkover detector in which the walkover detector is used to establish an overhead position directly above the boring tool 
using a locating signal transmitted at a single frequency. The measured overhead signal strength of the locating signal 
transmitted from the boring tool is then used in conjunction with the determined value of the skin depth to determine the 
depth of the boring tool below the overhead position on the surface of the ground such that the depth of the boring tool is 

30 established based at least in part on the skin depth. 

In another locating scenario, with the locating signal transmitted at a single frequency, the boring tool moves 
through the ground along an intended path while transmitting the locating signal and moves in an orientation which 
includes pitch. The boring tool includes pitch sensing means and the locating signal exhibits a field defined forward 
point at the surface of the ground with the boring tool at a particular point along the intended path. The field defined 
35 forward point being vertically above an inground forward point on the intended path through which the boring tool is 
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likely to pass. The boring tool is located by using a walkover detector to receive electromagnetic data which identifies 
the forward point. Signal strength of the locating signal is then measured at the forward point, as transmitted from the 
boring tool at the particular point, and the measured signal strength of the locating signal is used at the forward point in 
conjunction with the determined value of the skin depth and a sensed pitch value to determine the depth of the boring 
5 tool referenced to the particular point and to determine a forward distance on the intended path from the particular point 
at which the boring tool is located to the in-ground forward point. 

Alternatively, the field defined forward point may be located on or immediately above the surface of the ground 
and an overhead point may be identified on or immediately above the surface of the ground directly above the boring 
tool at the particular point. The forward distance is measured between the overhead point and the forward point as, 
10 defined at the surface of the ground. Using the forward distance, the determined value of skin depth and certain 
characteristics of the locating signal at the forward point, a skin depth corrected depth of the boring tool at the particular 
point is determined. 

In another alternative, the intended path of the boring tool in the region is configured such that the forward point 
is at a higher elevation on the surface of the ground than the particular point. The actual depth of the boring tool is then 
established at the particular point and a vertical elevation difference is measured between the particular point and the 
m forward point. Thereafter, the locating signal is sensed at the forward point while the boring tool is at the particular point 
W to determine an uncorrected depth of the boring tool which is subject to skin depth error. Using the measured vertical 
2 elevation difference, the actual depth of the boring tool at the particular point and the uncorrected depth of the boring 
hj tool measured from the forward point, a forward point skin depth correction factor is determined. During subsequent 
It) drilling operations the forward point skin depth correction factor is used in determining skin depth corrected depth with 
^ the boring tool at subsequent particular points. 

M BRIEF DESCRIPTION OF THE DRAWINGS 

Q The present invention may be understood by reference to the following detailed description taken in conjunction 

with the drawings briefly described below. 

25 FIGURE 1 is a diagrammatic illustration of a coordinate system for purposes of describing the well know 

magnetic dipole equations. 

FIGURE 2 is a diagrammatic elevational view of a horizontal drilling operation being performed in a region of 
ground using a portable walkover detector, shown here to illustrate the effects of skin depth on locating the boring tool. 

FIGURE 3 is a diagrammatic elevational view of another horizontal drilling operation being performed in a 
30 region of ground using a portable walkover detector, shown here to illustrate the effects of skin depth on locating the 
boring tool with the walkover detector at a forward locate point. 

FIGURE 4 is diagrammatic elevational view of still another horizontal drilling operation being performed in a 
region of ground using a locating/tracking system including fixed position above ground locating field detectors, shown 
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here to illustrate the effects of skin depth on locating the boring tool in a system using such fixed position above ground 
locating field detectors. 

FIGURE 5 is a plot of transverse and radial absorption parameters against the ratio of range to skin depth, 
shown here to illustrate the nature of these absorption parameters with increasing range to skin depth ratios. 

5 FIGURE 6 is a plot of absorption parameters versus the ratio of range to skin depth for rich agricultural earth 

extending to infinity in all directions, shown here to illustrate selected points on the absorption parameters for use in 
validating the multi-frequency approach of the present invention. 

FIGURE 7 is a block diagram illustrating a multi-frequency transmitter manufactured in accordance with the 
present invention and suitable for use in a boring tool. 

10 FIGURE 8 is a block diagram illustrating a multi-frequency receiver manufactured in accordance with the 

present invention and suitable for use in a portable walkover detector or fixed position detector. 

FIGURE 9 is a perspective view of the surface of the ground on which a calibration procedure is being 
performed for determination of skin depth in accordance with the present invention using a portable walkover detector. 

IH FIGURE 10 is a plot illustrating a two point calibration function F. 

55 FIGURE 1 1 is a diagrammatic elevational view of an overhead position determination setup, illustrating the 

IS determination of a skin depth corrected depth of the boring tool. 

Q FIGURE 12 is a plot representing the deviation from the cubic law obtained from a curve fit of Wait's 

r 1 theoretical results valid for a semi-infinite conductive region and zero transmitter pitch in the range 0 < D/8 < 3 . 

S FIGURE 13 is a diagrammatic elevational view of a forward locate position determination setup, shown here to 

130 illustrate determination of a skin depth corrected depth of the boring tool from the forward locate point. 

FIGURE 14 is a diagrammatic elevational view of a surface offset distance position determination shown here 
to illustrate a variation in which an above ground measurement between the forward locate point and an overhead point 
directly above the boring tool is used in determining the skin depth corrected depth of the boring tool. 

FIGURE 15 is a diagrammatic elevational view of an above ground elevation offset distance position 
25 determination setup, shown here to illustrate still another variation in which a forward locate point skin depth correction 
factor is developed for use in subsequent drilling operations. 

FIGURE 16 is a diagrammatic elevational view of a dual point common elevation position determination setup, 
shown here to illustrate the use of two identifiable points having a substantially common elevation. 



4 



DETAILED DESCRIPTION OF THE INVENTION 



Attention is immediately directed to Figure 2 which illustrates a boring system 10 operating in a region 12. It is 
noted that like reference numbers are used to refer to like components wherever possible throughout the various figures. 
The surface of the ground is indicated by the reference number 14. System 10 includes a boring tool 16 that is positioned 
on the end of a drill string 18 which is only partially shown. Boring tool 16 includes a dipole transmitter 20 having an 
antenna 22 that transmits a dipole locating field 24. The latter is received using a walkover portable locator/detector 30. 
Specific details regarding the implementation of system 10, as well as details regarding the implementation of other 
types of systems in accordance with the present invention will be given at appropriate points hereinafter. For the 
moment, discussions will be limited to more general details regarding skin effect as related to the operation of system 10 
in order to facilitate the reader's understanding. 

Still referring to Figure 2, locating field 24 , emitted by antenna 22 into a non-conducting medium (not shown) 
such as air is described mathematically by the cubic law, as mentioned above. However, drilling takes place in the earth 
such as, for example, in region 12 which is assumed to possess electrical conductivity characteristics. These 
characteristics result in significant deviations from the cubic law when locating field 24 varies sinusoidally in time at a 
frequency f. With conductivity of region 12 denoted as a, the penetration distance of locating field 24 as well as the 
shape of the magnetic field lines which make up the field depend on the frequency and conductivity parameters. 
Penetration distance is often called skin depth and is defined as: 



where (x denotes the permeability of the earth and 8 is the skin depth. Thus, skin depth decreases if conductivity, a, 
permeability, \i, or frequency, f, increase. Conversely, skin depth becomes infinite at zero frequency or conductivity, in 
which case the magnetic field is again described by the magnetic dipole relationship. The significance of the zero 
frequency characteristics will become apparent at an appropriate point below. Increasing conductivity and/or frequency 
serve to weaken the magnetic flux intensity recorded by locator 30 above ground. If skin effect is not accounted for, the 
boring tool can appear to be at a position farther from the detector than in actuality. In the present example, where 
measurements are being taken with a walk-over locator directly over the drill head (OH), the transmitter can appear 
deeper, at position A, where the boring tool is shown in phantom. 

Referring now to Figure 3, system 10 is shown once again with locator 30 at a different position. Specifically, 
the locator is shown at what is referred to as a forward negative locate point (FNLP) or, more simply, forward locate 
point (FLP) (see, for example, above referenced U.S. patent no. 5,337,002). An "X" is indicative of the configuration of 
the receiving antenna within locator 30. At the FLP, the flux lines of the locating field are characteristically vertically 
oriented. In the absence of skin depth compensation, when the locator is at the forward locate point, dipole transmitter 22 
can appear deeper and shifted farther away from the locator at position B where the boring tool is shown in phantom. 
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US patent application serial no. 08/835,834, filing date April 16, 1997, entitled Systems, Arrangements and 
Associated Methods for Tracking and/or Guiding an Underground Boring Tool is commonly assigned with the present 
application and is incorporated herein by reference. One example of a highly advantageous locating/guidance system 
conforming with the subject application is shown in Figure 4 and is generally indicated by the reference number 40. 
5 System 40 uses one or more above ground detectors 42 positionable at fixed locations within region 12 for reception of 
locating signal 24. With regard to the present application, skin depth produces an effect in system 40 which is similar to 
that described with reference to Figure 3. That is, boring tool 16 appears (shown in phantom) to be at position C at a 
deeper depth and shifted away from detectors 42. 

In some prior art systems, an above ground calibration procedure (not shown) is performed in an attempt to 

10 measure the signal strength of the dipole transmitter to be used in the boring tool with no consideration of the influence 
of skin effect. For example, the dipole transmitter and the locator are placed on the surface of the ground at a known 
separation and orientation. In this regard, it is recognized herein that skin depth has an effect on signal strength when 
such a calibration procedure is performed. Moreover, the skin effect in this above ground procedure varies from the skin 
effect encountered when the boring tool is beneath the surface, which further complicates provisions for skin depth 

15 compensation. The accuracy of skin depth obtained from an above ground calibration depends to a great extent on the 
J homogeneity of the soil. Skin depth will be accurate if the conductivity of the soil near ground surface, where the 
Cfi calibration has been performed, is similar to that of the soil above the transmitter. Specifically, in the calibration 

11 procedure for ranges up to about three times skin depth, skin effect causes an increase in recorded strength of the 
5 transverse component of flux intensity and hence the distance between transmitter and receiver appears smaller while an 
||p opposite trend can be observed for the radial component of the flux intensity. 

L Having generally described the influence of skin depth, attention is now directed to details concerning 

§pi provisions for effective compensation. The discussion immediately above, concerning an above ground calibration 
H : procedure, evidences that each component of the magnetic flux intensity is affected differently by earth conductivity and 
£ dipole frequency. An exact solution, available for the components of magnetic flux intensity of a magnetic dipole 
US immersed in homogenous earth (i.e., having a uniform conductivity) of infinite extent is given in US patent no. 
4,710,708, issued to Rorden et al. Rorden, however, relies on the solution only to show that skin effect can be ignored 
when the range of interest is significantly less than the skin depth. Accordingly, Rorden uses a locating frequency that is 
low enough (generally 1-100 Hz) to produce a sufficiently high skin depth in equation 5 above to accomplish this 
objective. The present invention, however, considers the use of such low frequencies as unacceptable because common 
30 signal detection hardware and sensors such as coil loops are more sensitive at higher frequencies. Additionally, state of 
the art systems such as, for example, systems 10 and 40 described above contemplate the use of the locating signal as a 
carrier for the purpose of transmitting data to above ground locations wherein the data are encoded upon the locating 
signal Carrier frequencies in the range of 1- 100 Hz limit data transmission capabilities optimistically to rates in the 
range of only 0.5 to 50 Hz according to the Nyquist criteria. 

35 Referring to Figure 5, transverse and radial absorption parameters are shown plotted against the ratio of range, r, 

to skin depth, 8, based on the Rorden expressions for these values. Exponential decay is shown for comparison. It is of 
interest here to note that Rorden fails to provide a plot such as that of Figure 5 and, even though Rorden lists the 
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equations, he fails to observe that the transverse and radial components are each affected in a different way as r/8 
increases. Both components, however, approach exponential decay asymptotically for large values of r/5. At the same 
time, it is important to understand that the usefulness of the Rorden teachings is inherently limited, in a practical sense, 
because the presence of a ground surface is not modeled. As will be seen, the present invention resolves the difficulties 
5 in using locating signals having frequencies sufficiently high for providing ease of detection and adequate data transfer 
rates using the locating signal as a carrier while providing effective and highly advantageous compensation for skin 
effect, even though the locating frequencies used are high enough to encounter levels of skin effect that should not be 
ignored. 

In one implementation of the present invention, which is applicable to essentially any underground boring 
10 system including systems 10 and 40 described above, the recognition must be emphasized that the locating field is 
governed by the cubic law of a magnetic dipole if the signal frequency goes to zero, because the skin depth goes to 
infinity regardless of the conductivity of the earth. At first blush, this recognition may seem of little importance to one 
of ordinary skill in the art since, as a practical matter, static magnetic fields (i.e., at zero frequency) are useless in the 
present application. That is, only a time varying field is readily measurable with sensitive detectors such as coils, 
flf However, the present invention overcomes the seeming uselessness of attempting a solution at zero frequency by 
*0 . providing a highly advantageous multi-frequency approach which allows the formulation of a zero frequency solution. 
Yi Moreover, the multi-frequency approach taken by the present invention may readily be implemented using existing 
.5 technologies, as will be seen immediately hereinafter. 

yj The multi-frequency approach of the method of the present invention requires transmission and reception of the 

S lb locating field using at least two different frequencies. While specific implementations to be described rely on the use of 
O four different frequencies, it is to be understood that any number of frequencies of two or more may be employed. To 
111 some extent, it is considered that accuracy may be enhanced, however, when more than two frequencies are used. Details 
^~ regarding frequency selection will be provided at an appropriate point below. In the four frequency implementation 
O under discussion, measurement of the locating field is made at one or more above ground locations using either one or 
fe more detectors configured for use at fixed locations and/or using a portable walkover locator. Receivers in either a 
portable walkover locator or in fixed position above ground detectors may be configured in essentially the same manner, 
in accordance with these teachings. During operation of any system utilizing the multi-frequency approach of the 
present invention, measurements are made at the above ground locations corresponding to each of the selected 
frequencies. Thereafter, these measurements are utilized in a highly advantageous way which serves to extrapolate a zero 
30 frequency measurement. For example, it may be assumed for any particular above ground location that the component 
of the magnetic field intensity of the locating signal measured at the i-th frequency £ depends on skin depth according to: 



(6) S,=S 0 F 
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where S 0 is a constant which corresponds to the intensity of the locating field at zero frequency, 6 ; is the skin depth at 
each of the selected frequencies, D is the depth or some characteristic length scale of the boring tool and F is a function 
to be determined. Thus, the objective is to establish the value of S 0 based on the values S,. To that end, for each of the 
selected frequencies, an interpolation polynomial or any other suitable mathematical function including an exact 
solution, if obtainable, may be used to provide a curve fit to the measured data for each of the frequencies. As an 
example, a cubic polynomial can be used to approximate the function, F, at the four required frequencies. Introducing 
the definition of skin depth, the magnetic field intensities can be written as: 

(7) S^S. + af^ + bf. + cf^ 

(8) S 2 = S 0 + af 2 0 5 +bf 2 + cf 2 15 

(9) S 3 = S 0 + af 3 0 5 + bf 3 +cf 3 

(10) S 4 = S 0 + af A 05 +bf A + cf A 



1 5 

1 5 



Equations 7-10 are a set of linear equations for the unknown coefficients S 0 , a, b and c that can be solved 

jj) employing standard solution methods. It should be noted that this approach is very efficient numerically, requiring a 

55 small matrix to be inverted with coefficients depending on the chosen frequencies. Once a value for S 0 is obtained, the 

W position of the boring tool can be determined using the well known cubic equations 1-4 above. Remarkably, there is no 

£ need to determine the skin depth values 5,. Since the selected frequencies are chosen prior to initiation of drilling, the 

|Il inversion of this matrix need only be performed once. Other formulations of signal strength at each frequency, S 19 may 

T5 be used. For example, one possible formulation may be based on the exact solution of Rorden in which exponential 

Q decay is observed for large values of r/S, the equations may be expressed in the form: 



(11) S r = S 0 e *<G^ 

(12) S, = e-«"(S 0 +afr+bf,) 



where G is a function to be determined and all other values are described above. It should be noted that this 
approximation also requires measurements at four frequencies but the four resulting equations for unknown coefficients 
20 S 0; a, b, c are nonlinear. Hence, the solution method is somewhat less efficient than one based on polynomial 
approximations, but remains applicable over a wider range of skin depth. Once again, there is no need to determine the 
values Sj. 

Referring to Figure 6, an important feature of the multiple frequency approach of the method of the present 
invention resides in the fact that it does not require knowledge of earth properties or ground surface geometry. The 
25 components of the magnetic field measured at nonzero frequencies contain property and geometry effects and pass them 
on to extrapolated zero frequency values. Figure 6 demonstrates the validity of this approach for a simplified case, rich 
agricultural earth (i.e., meaning higher conductivity soil) extending to infinity in all directions. In this example, any 
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variations of earth conductivity and the effect of ground surface and air on the magnetic field are neglected. The 
transverse and radial absorption parameters, representing the deviation of the magnetic field from the cubic law of 
magnetic dipoles, are plotted vertically against the ratio of range to skin depth as dashed lines indicated by the reference 
numbers 50 and 52, respectively. Points 50a-d and 52a-d have been selected on each absorption parameter curve 
5 corresponding to transmitter frequencies of 2, 8, 14 and 20 kHz, respectively. As seen, extrapolation of the curve using 
the polynomial formulation of equations 7-10 provides magnetic field data that are within 2% of known exact values for 
the absorption parameters at zero frequency (i.e., the known value of 1.0 for both of the absorption parameters at zero 
frequency). In the method of the present invention, the extrapolation is performed for magnetic intensity, however, the 
present example serves to illustrate the validity of this approach even though the absorption parameters were 
10 extrapolated since magnetic intensity is the product of the cubic law and absorption parameter. Even though the 
discussion dealt with distributed conductivity of soil, the multiple frequency approach will work for other field 
distortions due to conductivity including but not limited to buried electrical conductors, pipes, plates and rebar. 

With regard to selection of frequencies at which the locating signal is to be transmitted, it is noted that an 
unlimited number of different frequency combinations may be employed. Since an extrapolation to zero frequency is 
being performed, however, it is considered that the frequencies should be as low as possible while still providing for 

5 adequate detection and transmission of data, using the locating signal as a carrier. For example, frequencies in the range 

fi of 2-40 kHz are considered as acceptable. 

% Referring now to Figure 7, having described the multi-frequency approach of the method of the present 

W invention, descriptions will now be provided of components appropriate for use in systems which utilize the approach. 
%) Figure 7 illustrates a multi-frequency transmitter manufactured in accordance with the present invention and generally 
f== indicated by the reference number 100. Transmitter 100 includes a sensor/conditioning section 102 and a carrier 
m generation/antenna drive section 104. Transmitter 100 is generally configured for use in a boring tool, in certain 
instances, the transmitter may be used in above ground applications such as, for example, in a calibration and/or system 
Q test unit (not shown). 

25 Still referring to Figure 7, sensor/conditioning section 102 includes a suitable group of sensors in this instance 

comprising, for example, a three-axis gyro 106, a three-axis magnetometer 108, a three-axis accelerometer 110, a 
roll/pitch sensor 112, a temperature sensor 114 and a battery sensing section 116. Physical parameters at the outputs of 
magnetometer 108, accelerometer 110 and roll/pitch sensor 112, as well as the transmitter battery conditions using 
battery sensing section 1 16 and temperature using temperature sensor 1 14, are provided to a multiplexer 118 which then 

30 transfers all of these signals in multiplexed form to an analog to digital converter 120. The latter digitizes and converts 
the multiplexed signals into digital format, for example, at either an 8-bit or 12-bit resolution, depending on accuracy 
requirements. Thereafter, a microprocessor 122 processes all of the parameters provided from the analog to digital 
converter and converts the parameters into information relating to the in-ground transmitter coordinates or relating to 
down-hole conditions. For example, during this parameter processing, the microprocessor may perform linearization and 

35 temperature compensation on the output of roll/pitch sensor 1 12 in order to calculate an absolute pitch position of the 
boring tool. The pitch output may further be compensated based on the determination of a particular roll position. 
Linearization and compensation coefficients are generated during factory calibration and stored in a calibration constant 
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section 124 which comprises a non- volatile memory area within transmitter 100. It is noted that programming of 
microprocessor 122 is considered to be within the ability of one having ordinary skill in the art in view of this overall 
disclosure. 

Once the calculations are complete relating to all of the signals from sensors in sensor/conditioning section 102, 
5 the results are transmitted to one or more above ground locations. To that end, carrier generation/antenna drive section 
104 includes an oscillator 126 which provides a clock signal to microprocessor 122 and to a divide by N counter 128. 
The latter receives a frequency select input from microprocessor 122 on an f select line 130 such that the divide by N 
counter may selectively generate any one of a wide range of carrier frequencies. The carrier frequency may be selected 
by microprocessor 122 under software control in a number of different ways. In one configuration, frequency selection 
1 0 can be performed at the beginning of the drilling operation, for example, after monitoring of background noise levels at 
various frequencies such that noisy frequencies may be avoided. Such noisy frequencies may be attributed, for instance, 
to traffic loops, invisible dog fences, cable TV, and power lines in a particular region. In another configuration, the 
transmitter may change frequencies on the fly, after the drilling has started. On the fly frequency change can be initiated 
either by the microprocessor using a pre-determined algorithm, or by the request of the drilling operator, for example, 
JL5 using a signal transmitted by telemetry from the surface to the boring tool or transmitted through the drill string using an 
isolated electrical conductor or based on possibilities such as, for example, boring tool roll orientation sequence and roll 
03 rate. Particularly advantageous arrangements for automatically forming an isolated electrically conductive path between 
1~ the drill rig and an in-ground device such as a boring tool to provide power and signal paths are disclosed in co-pending 

p U.S. patent application serial no. 9 filing date (attorney docket no. DCI- 

P0 1 8), which is incorporated herein by reference. 

The selected carrier frequency is then passed to a modulation section 132 which is configured for modulating 
01 data from sensor/conditioning section 102 onto the selected carrier frequency. Modulation section 132 receives the data 
^; on a data line 134 from microprocessor 122 and also receives a modulation selection signal on a mod select line 136 
S connected with the microprocessor. The modulation selection signal may select, for example, phase modulation or 
|35 amplitude modulation, or combinations of both. The modulation scheme may be programmed either before or during the 
drill, much in the same manner as in the case of the carrier frequency, described above. 

With continuing reference to Figure 7, data modulation section 132 passes a modulated carrier signal to a driver 
section 140. The driver section receives a power selection input on a power select line 142 from microprocessor 122. In 
this way, the output of transmitter 100 may be tailored to drilling conditions, for example, to conserve battery power in a 

30 shallow drill run or to increase transmitter output at longer ranges and/or drilling depths or even to stop transmission 
altogether during idle periods of a drilling operation. Control of the power as well as other functions can be achieved 
using procedures such as have been described for frequency control. The transmitter can send encoded data (as is done 
for roll, pitch and other parameters) to allow the receiver to adjust its calibration for the new signal strength. This will 
allow the operator to continue monitoring depth or range without the need to recalibrate while drilling. An antenna drive 

35 signal is produced on an antenna line 144 which is coupled to an antenna 146 which generally comprises a dipole 
antenna for emanation of locating signal 24. 
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In accordance with the multi-frequency approach of the present invention, described above, transmitter 100 may 
alternately transmit one of four selected carrier frequencies from the boring tool. The carrier frequencies may alternate at 
any suitable rate such as, for example, 10 Hz and may be selected in accordance with considerations described 
previously. It should be appreciated that transmitter 100 is configured for flexibility in carrying out the method of the 
5 present invention. That is, fewer or more than four carrier frequencies may readily be transmitted either individually or 
simultaneously. 

Turning now to Figure 8, a multi-frequency receiver manufactured in accordance with the present invention is 
generally indicated by the reference number 150. Locating signal 60 is received by an antenna arrangement 152 which 
may include three orthogonally arranged x, y and z antennas indicted by reference numbers 154, 156 and 158, 

10 respectively. The locating signal received by each antenna is first amplified by respective very low noise pre-amplifiers 
160a-c prior to processing. This pre-amplification maintains the received signal-to-noise ratio while making any noise 
introduced by subsequent circuitry relatively negligible. The amplified antenna signals (not shown) are then fed to 
mixers 162a-c to be translated down to a lower intermediate frequency (IF). In this manner, any locating signal 
frequency within a suitable design range such as, for example, 2-40 kHz can be received through appropriate 

15 adjustment of the frequencies of the mixers such that outputs of the mixers fall within a selected IF band. 

m Still referring to Figure 8, the mixers are followed immediately by narrow-band, band-pass filters 164a-c which 

V) essentially pass only the locating signals at the translated IF frequency. The IF frequency can be either the sum or the 

if* 

X. difference of the carrier frequency and the mixer frequency. The filtered x, y and z signals are further amplified by 
k! programmable-gain amplifiers (PGA's) 166a-c before being received by an analog-to-digital converter 168. PGA's 166 
lb provide over 96 dB of dynamic range and are each directly controlled by a digital signal processor (DSP) 170. One 
J- suitable DSP is the ADSP2185L, a sixteen bit fixed point DSP manufactured by Analog Devices, Inc. Analog-to-digital 
01 converter (ADC) 168 digitizes the received signals, at a rate controlled by a direct digital synthesizer (DDS) 172 which 
is, in turn controlled by DSP 170. The DSP and DDS receive an oscillator signal from an oscillator 174. Using the 
oscillator signal and based on control from the DSP, the DDS generates a local oscillator frequency (LO) for mixers 162. 
Ii5 The ADC digitizes the received signals at the rate determined by the DDS and converts the signals to a binary number 
two's complement format. The conversion rate is either four times the IF frequency, if a quadrature sampling scheme is 
used, or may be significantly less than the IF frequency, if an under-sampling scheme is chosen. The resolution of the 
ADC may be 12-bit to 16-bit. All axes are simultaneously sampled in order to maintain relative phase. 

Continuing to describe receiver 150, digital signal processor DSP 170, controls all operations of the receiver 
30 including mixing frequency, PGA gain, and a selected signal processing algorithm. In the case of a quadrature sampling 
scheme, the DSP samples the received signals at four times their IF (translated) frequency and then multiples the 
received signals it by a separate Sine and Cosine sequence to obtain in-phase and quadrature-phase components. This 
process converts the received signal from its IF frequency down to a base-band frequency that contains modulated data, 
if present, while, at the same time, breaks down the signal into its in-phase (I) and quadrature-phase (Q) components. 
35 The I and Q components are each passed through a simple low-pass filter (not shown) to remove everything but the 
modulated data. The filtered outputs are then used to obtain the original data as well as further processing to recover 
signal magnitude and sign information. Additionally, the outputs are also used, along with a modified phase-lock-loop 



11 



technique known as Costas loop, for controlling the DDS frequency (which controls the mixer frequency and ADC 
sample rate) and the PGA gain settings. The exact algorithm varies depending on the modulation scheme used but may 
be developed by one having ordinary skill in the art in view of this overall disclosure. 

If under sampling (not shown) is used, the DSP would sample the received signals at a rate much lower than the 
5 IF frequency. The digitized data is then processed using a matched filter to obtain data, magnitude, and sign information 
as well as for PGA gain control. Irrespective of sampling, the DSP implementation is considered to be highly 
advantageous, resulting in a very flexible and adaptive multi-frequency receiver. Many modifications (not shown) are 
possible in view of this disclosure for purposes of performance improvement. For example, mixers 162 can be eliminated 
by replacing narrow-band band-pass filters 164 with broadband band-pass filters and using ADC 168 to perform 
10 quadrature sampling and direct-to-base-band conversion (digital mixing) in a single operation. In order to receive the 
locating signal at different frequencies, DSP 170 may either sample the data at different rates or sample everything at a 
single, fixed rate and then perform rate conversion in software using decimation and interpolation techniques (known as 
digital re-sampling). As mentioned previously, it should be appreciated that receiver 150 may readily be incorporated 
into either a portable walkover locator or into detectors designed for use at fixed positions with a drilling region. 

4tf5 Having described a highly advantageous multi-frequency approach for use in skin depth compensation, skin 

m depth compensation techniques using a single frequency locating signal will now be described with regard to a number 
111 of different exemplary scenarios. It is to be understood that existing systems using portable locators may readily be 
E adapted in conformity with these teachings or, alternatively, new systems using either a portable locator and/or one or 
hj more locating field detectors designed for positioning at fixed locations within a drilling region are also readily adaptable 
^0 in view of these teachings. 

r2 Referring now to Figure 9, the techniques disclosed herein for use with single frequency locating signal 

transmission initially rely on a determination of the skin depth in the drilling region. As mentioned previously, some 
4^ prior art systems utilize an above ground calibration procedure in an attempt to relate the signal strength of a dipole 
T\ transmitter to distance without skin depth compensation. The present invention introduces a technique for performing an 
25 above ground calibration procedure which not only provides dipole signal strength, but also yields a value for skin depth 
in the drilling region which may be used in subsequent position determination techniques accounting for skin depth. 
Figure 9 illustrates a calibration procedure being performed using a portable walkover detector 260 in a region 262. The 
calibration procedure is performed on the surface of the ground which is assumed to be planar for purposes of simplicity, 
having x and y coordinate axes defined as shown. A dipole transmitter 264 is diagrammatically illustrated and is oriented 
30 along the y axis while being centered upon the x axis. Preferably, the transmitter should be the transmitter which is to be 
used during subsequent drilling operations in a drilling configuration such as housed in the drill head (not shown) placed 
on the surface of the ground. Alternatively, the transmitter itself may be positioned on the ground, but it must be 
remembered that measurements are likely to be affected by any housing later positioned around the transmitter. 

Still referring to Figure 9, the calibration procedure is performed with the walkover locator at two offset 
35 positions along the x axis indicated by the reference numbers 266a, corresponding to an offset distance of xl, and 266b, 
corresponding to an offset distance of x2. Transmitter 264 transmits locating signal 24 at a single frequency. It is to be 
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understood that the calibration procedure may just as readily be performed using a detector which is intended for location 
at a fixed position within the drilling region following the calibration procedure. In this regard, irrespective of the 
specific form of the detector instrument to be used in the calibration procedure, the instrument should be positioned such 
that its locating field sensor arrangement is on the x axis. In the instance of a walkover locator having an antenna 
5 configuration as described, for example, in above referenced U.S. patent no. 5,337,002, which is commonly assigned 
with the present application, the plane of the antenna arrangement should be aligned parallel to the y axis of the 
transmitter. In the instance of a locating field detector including three orthogonal receiving axes, such as described in 
above incorporated US application serial number 08/835,834, the detector arrangement is somewhat arbitrary since 
signals measured along the three axes can be transformed mathematically into any desired directions. 

10 Using the configuration shown in Figure 9, the calibration procedure is performed by measuring the 

components of the magnetic flux intensity B yl and B y2 at positions Xi and x 2 , respectively. Geophysical theory provides 
an equation for the calculation of dipole strength and skin depth that has the general form: 

(13) B y = B y (x,y,w,S) 

^ where B y is a measured intensity, x and y are the coordinates of the locator/detector, m is the signal strength of the dipole 
JE5 transmitter and 5 is the skin depth. At this time, a preferred method is based on the theory of Wait et al (Journal of 
Geophysical Research, Vol. 58, No. 2) which is valid for zero transmitter pitch, level ground surface and homogeneous 
= f. soil conditions, as are present in Figure 7. Wait solves Maxwell's equations with boundary conditions at the ground 
? surface correctly satisfied. Since the calibration procedure provides two values of magnetic flux intensity and the 

fS distances xi and x 2? two nonlinear equations for calculating m and 5 are obtained: 

_£ (14) B y] =B y (x,Am,S) 

E (15) B y2 = B,(x 29 0,m 9 S) 

20 

The deviation of B y from the cubic law is approximated in the range 0 <= x/8 < 3 by: 



(16) 



(17) 



*§-.♦»§♦«(* 



+ d\ 



where the function F is shown in Figure 10. The unknown coefficients b, c and d can be obtained from this graph using 
standard numerical techniques. 
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Using equations 16 and 17, the following equations can be obtained which can be solved for m and 8: 



(18) 
(19) 



r - m 



R - m 

X 2 



The solution is obtained in 2 steps. First, the following variables are defined after introduction in equations 18 

and 19: 

(20) g, = B vl x, 3 

(21) g 2 = £ v2 x 2 3 

(22) * = i 



Subtracting equation 19 from equation 18 provides equation 23 for 

(23) ( gl x 2 - g 2 x*)de 3 + (g } x 2 -g 2 x x 2 )cs 2 +(g,x 2 -g 2 x l )be + g ] - g 2 = 0 



Equation 23 can be solved employing a standard method such as Newton's to yield 5. Thereafter, dipole 
strength, m, follows directly from: 



(24) 



m ■ 



B y\ X \ 



S 



+d 



Thus, dipole signal strength, m, and skin depth, 8, are established for use in subsequent position determinations. 

Turning now to Figure 1 1, in a first scenario, an overhead position determination setup (hereinafter OH setup) is 
generally referred to by the reference number 270 with a detector (not shown) at a position 272 located directly overhead 
(hereinafter OH setup) of transmitter 264 transmitting locating signal 24. At this location, the flux lines of the magnetic 
locating field are characteristically horizontal substantially over the transmitter. The detector measures the horizontal 
component of the magnetic flux intensity B x . With m and 8 known from the foregoing above ground calibration 
transmitter depth and with a measured value of intensity, B xD , from the detector, D, is determined from a single equation 
written symbolically as; 
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(25) B xD = B r \m,D,^ 



The exact form of equation 25 can either be obtained from geophysical theory or from dimensional analysis. 
Applying the latter (e. g. } P. W. Bridgman, Dimensional Analysis, 1931) six variables are identified governing the 
physics of OH setup depth measurement. The variables include B s , ju, c, f, D and m which have been defined previously. 
Furthermore, four fundamental units including length, time, volt, and ampere characterize the problem. Hence according 
to the 7i -theorem of dimensional analysis six minus four or two non-dimensional groups describe the OH setup 
measurements mathematically. The two non-dimensional groups include 



(26) 



m 



(27) fi cx/D 2 

The second group given by equation 27 can be simplified to D/8 using the definition of skin depth from 
equation 5. Hence B x must be of the following general form: 

Here, the function G represents the deviation from the cubic law obtained from a curve fit of Wait's theoretical 
results valid for a semi-infinite conductive and zero transmitter pitch in the range 0 < D/5 < 3 medium, shown in Figure 
12. Since this equation is nonlinear for depth D, an iterative procedure must be formulated. As one example: 



(m 


G 




U 







(29) D = 



Function iteration/successive approximation is performed beginning with an initial guess for D, e. g. the value 
corresponding to infinite skin depth. In successive approximations, the procedure inserts the last available value for D on 
the right hand side of this equation thereby calculating a new, more accurate value. This process is repeated until 
changes between successive values of D are reduced to a specified tolerance. 

The analysis outlined immediately above provides the correct functional relation between variables governing 
OH setup depth measurement which can be written as: 
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The unknown coefficients c } and d % must be obtained from another source, for example, Wait's theory or a 
physical experiment conducted in different soil conditions and at various depths. Another method for obtaining these 
coefficients relies entirely on numerical modeling solving Maxwell's equations and pertinent boundary conditions. 
5 Computer codes are commercially available to aid in this task such as, for example, software by Infolytica Corporation in 
Montreal, Canada. 

Attention is now directed to Figure 13 which illustrates a second scenario representing a Forward Locate Point 
position determination setup generally indicated by the reference numeral 290. It should be mentioned that even though 
the present discussions are made with reference to the forward locate point, these concepts are equally applicable to the 
10 rear locate point In this instance, detector 30 measures the magnitude of the magnetic flux intensity B (not shown) as a 
result of locating signal 24 transmitted from transmitter 264 within drill head 274. In addition, transmitter pitch y is 
measured. In order to calculate transmitter depth from these measured quantities the following equations must be solved: 



(31) B { = B 4 \t,rj, s 



(32) B n = B n \4,r,, s 

B £ 

(33) Um r = -± 

(34) D = r sin(2 + y) 

(35) B y = VV + V 

(36) tan A = ^~ 

V 

(37) r 2 = ^ + 7 2 



The definitions of the geometric variables D, r, r|, X are given in Figure 13. It is noted that transmitter depth is 
15 the distance from the ground surface to the transmitter. Since the locator antennas measure signals above ground, the 
distance between the antennas and ground must be subtracted from the computed depth. These variables and the two 
components of the magnetic field intensity and B n make up a total of 7 unknowns that can be obtained from the listed 
seven equations using standard numerical methods. An example of a convenient solution method is to rewrite the 
equations in terms of polar coordinates r, X using 

20 (38) | = rcos/t 

(39) rj = rsin/L 

This transformation eliminates two equations. The remaining equations then read: 
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(40) 



(41) 



B < = B { [ r >*>j) 



3 cos 2 X - 1 



3 sin X cosX 



r 



r 




(42) 



tan y = — — 



(43) 
(44) 



D = f sin(/l + /) 



Note that the equations for the components of magnetic field intensity express the cubic law of a magnetic 
dipole multiplied by a function H x or H 2 that accounts for the effect of skin depth. The latter is known from an above 
ground calibration, as described above. Details of these functions can be derived employing either geophysical theory or 
dimensional analysis. Further, note that equation 44 for transmitter depth is uncoupled from the other equations allowing 
independent solution for the position coordinates r and X of the FLP based on the following nonlinear equations: 



where the variables have been defined above. These equations can be solved employing any of the standard solution 
methods for sets of nonlinear equations such as, for example, Newton's and function iteration. 



generally indicated by the reference number 320. This setup technique uses the FLP and OH points indicated by the 
reference numbers 322 and 324, respectively. In this regard, it is noted that the locations of these points are affected as a 
result of skin depth. At forward locate point 322, the horizontal component, B x , of the magnetic flux intensity vanishes. 
This fact can be used to derive a formula for transmitter depth D as a function of a horizontal distance, Ax G , at the 
surface of the ground between OH point 324 and FLP 322. In applications where skin effect can be neglected, a simple 
equation for depth, D, can be derived from the cubic law for a dipole field: 



In order to account for skin effect relying on B x = 0, a different form of the equation is used which is written in 
symbolic notation as: 




Turning now to Figure 14, in a third scenario, a surface offset distance position determination setup is 



(47) 



D = j2Ax G 



(48) 




Here, skin depth 5 is obtained from an above ground two-point calibration as described earlier and Ax G can be 
measured easily using available standard distance measurement methods. Details of this equation can also be derived 
from geophysical theory, e. g., the aforementioned work published by Wait. In general, an explicit formula for depth 
cannot be derived from this equation since it will most likely be nonlinear in D, therefore, the expression must be solved 
numerically employing a suitable standard solution method such as Newton's or function iteration. 

Attention is now directed to Figure 15 which illustrates a fourth, above ground elevation offset distance position 
determination setup generally indicated by the reference number 360. This setup technique is useful in conjunction with 
an OH measurement of drill head depth, D 0H5 , which accounts for skin depth such as, for example, described above in 
the OH position determination setup associated with Figure 1 1 . The elevation offset technique requires a measurement of 
elevation change, Az G , between an OH point 362 and a FLP 364 with transmitter 264 at one position 366. Generally, this 
measurement will be performed once early in a drilling operation. With the drill head at position 366, FLP-depth, D F lps, 
( wherein 5 indicates compensation for skin depth) can be calculated from the over head depth, D 0 hs ? and the measured 
elevation change using: 

(49) D FLPS = D 0Hg + te G 



Still referring to Figure 15, the present technique is especially useful for a walk-over locator that is able to 
accurately measure over-the-head depth accounting for skin effect, but is not configured for skin depth compensation 
from the FLP. Since depth measured with such a locator at a forward locate point does not include skin effect, the ratio 
D FLP5 / D FLP obtained at transmitter position 366 can be employed to correct subsequent FLP depth measurements using 
the formula; 



(50) 



D 



r D OHS + Az G y 



FLP8, 



D 



FLP 



D 



FLPk 



Here, the subscripts 1 and k denote the first and k-th locating positions respectively. It should be noted that this 
correction is linear and hence can only be expected to give accurate results for small depth variations and homogeneous 
soil properties. 

Referring to Figure 16, a dual point common elevation position determination setup is generally indicated by 
the reference number 380. Accordingly, magnetic measurements are made at any two identifiable points having a 
common elevation. These measurements are combined to obtain transmitter depth in the presence of skin effects. It is 
assumed that ground elevation changes only moderately between these two points so that the same locator elevation can 
be maintained by simply raising or lowering the locator/detector unit relative to the surface of the ground. This approach 
does not require an above ground calibration since skin depth is determined from magnetic field data together with 
transmitter depth. Two such identifiable points are the OH and the FLP points, as indicated. At the OH location there is 
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only a horizontal magnetic flux B x and at the forward locate point there is only a vertical flux component B y . Even 
though transmitter 264 is shown in a level orientation in this figure, its pitch is not required to be zero. The governing 
equations have already been described above for techniques relying on separate OH and FLP measurements. 
Specifically, equations including (28) and (31) to (37) may be used for the solution of transmitter depth D, skin depth 5 

5 and the variables r 9 £ 7/, A, B^ defined above. An example of a practical solution method can be derived by 

modifying the approach given previously for solving the set of equations (31) to (37). There the method requires the 
simultaneous solution of the two nonlinear equations (42) and (43) for the geometric parameters X and r. Here D, 
8, X and r are obtained by solving the 4 nonlinear equations (28), (44), (45) and (46) simultaneously employing any 
standard numerical solution methods such as Newton's and function iteration. 

10 In that skin depth compensation arrangements and associated methods disclosed herein may be provided in a 

variety of different configurations and modified in an unlimited number of different ways, it should be understood that 
the present invention may be embodied in many other specific forms without departing from the spirit of scope of the 
invention. Therefore, the present examples and methods are to be considered as illustrative and not restrictive, and the 
invention is not to be limited to the details given herein, but may be modified within the scope of the appended claims. 
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